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Abstract 

The dynamical behavior of the magnetism of diluted magnetic semiconductors (DMS) has been 
investigated by means of atomistic spin dynamics simulations. The conclusions drawn from the 
study are argued to be general for DMS systems in the low concentration limit, although all 
simulations are done for 5 % Mn-doped GaAs with various concentrations of As antisite defects. The 
magnetization curve, M(T), and the Curie temperature Tq have been calculated, and are found to 
be in good correspondence to results from Monte Carlo simulations and experiments. Furthermore, 
equilibrium and non-equilibrium behavior of the magnetic pair correlation function have been 
extracted. The dynamics of DMS systems reveals a substantial short ranged magnetic order even at 
temperatures at or above the ordering temperature, with a non-vanishing pair correlation function 
extending up to several atomic shells. For the high As antisite concentrations the simulations show 
a short ranged anti-ferromagnetic coupling, and a weakened long ranged ferromagnetic coupling. 
For sufficiently large concentrations we do not observe any long ranged ferromagnetic correlation. 
A typical dynamical response shows that starting from a random orientation of moments, the 
spin-correlation develops very fast (~ Ips) extending up to 15 atomic shells. Above ~ 10 ps 
in the simulations, the pair correlation is observed to extend over some 40 atomic shells. The 
autocorrelation function has been calculated and compared with ferromagnets like bcc Fe and 
spin-glass materials. Wc find no evidence in our simulations for a spin-glass behaviour, for any 
concentration of As antisites. Instead the magnetic response is better described as slow dynamics, 
at least when compared to that of a regular ferromagnet like bcc Fe. 

PACS numbers: 75.25.-|-z,85.75.-d 
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I. INTRODUCTION 



The interest in dilute magnetic semiconductors (DMS) has been enormous during the 
last decade.-i^i^jii^iSiiliSiiSiiSiiiii^iMiM^ii^iiliiSiiiE Numerous theoretical and experimental inves- 
tigations focusing on their properties have been performed, motivated not least by the quest 
for a DMS system with ferromagnetic properties at and above room temperature. One of 
the more frequently studied materials is the III-V semiconductor gallium arsenide, doped 
with manganese (Mn-GaAs), and it now stands clear that an ordering temperature of 170 
K- to 180 can be obtained. Among the achievements on the theoretical side is the 
success of calculations based on density functional theory, where quite generally for DMS 
systems one should mention the following: first principles calculations result in 'reasonable' 
electronic band structures^ as compared to experiments, the correct values of the mag- 
netic moments are obtained^iii^ii^, calculations of interatomic exchange interactionsi^ii^ii^ 
result in ordering temperatures in agreement with observations, and finally the recognition 
of percolatioii^'^'^ as an important mechanism to govern the presence or absence of long 
ranged magnetic ordering, so that mean field theories should be used with great caution. 

Whereas most authors have concentrated on the static properties in thermal equilibrium, 
the scope of this Article is atomistic simulations of the dynamic properties of diluted mag- 
netic semiconductors. We have chosen to use Mn doped GaAs, with and without defects 
such as As antisites, as a representative of this group of materials and in our manuscript 
we argue that our findings apply to DMS in a quite general way. The main reason behind 
this is that both III-V and II- VI based DMS have very similar behavior as regards the 
exchange interaction, with an exponentially decaying behavior with distance between mag- 
netic impurities.'i^ In addition, for both sets of systems a varying degree of ferromagnetic 
and antiferromagnetic interactions can be obtained, depending on vacancy concentration of 
the host lattice. For instance, Mn doped ZnO has antiferromagnetic interactions which can 
be turned to ferromagnetic interactions in the presence of Zn vacancies.- In a similar way, 
it is known that Mn doped GaAs has primarily ferromagnetic interactions between the Mn 
atoms in the absence of defects, but an increasing amount of As antisites and Mn intersti- 
tials introduce antiferromagnetic interactions between the Mn atoms.— The main difference 
between the III-V and II- VI materials is that the solubility of transition metals is larger 
for the latter, but in the low concentration limit (2-10 % of transition metals doped into 



3 



the host lattice) these two sets of materials should have a similar magnetic behavior and 
dynamic response. 

Many DMS materials are in experiments analyzed as changing from a ferromagnetic 
behavior to a regime where spin-glass behavior is thought to occur. The strongest evidence 
for spin-glass properties can be found among the II- VI DMS systems, like Zni_a.Mna;Te^ 
and Zni_^.Mn^.In2Te42^'23. Among the III-V DMS materials spin-glass properties have been 
suggested for GaMnN— and Te doped GaMnAs^^. Often this transition is found as a function 
of concentration of magnetic dopants, but the concentration of other defects have also been 
analyzed to drive a spin-glass behavior. A good example of this is the system of interest here, 
Mn doped GaAs, where defects like As antisites, have been argued to produce a spin-glass 
like behavior^-. Spin-dynamics simulations, as presented here, are a natural way to analyse 
this behaviour and part of this study is devoted to analysing spin-glass behaviour in Mn 
doped GaAs. 

In the present paper we give a full account of a spin-dynamics simulation of a DMS 
material, including an analysis of short versus long ranged correlation and typical time 
scales of the magnetic response. We have considered Mn^^Gai^^'As with a Mn concentration 
of X = 5 % and a varying degree y = 0.00 — 2.00 % of As antisite concentration (As atoms 
which are positioned on the Ga sub-lattice). As this paper reports atomistic spin-dynamics 
simulations of DMS materials we note that one previous spin-dynamics simulation has been 
pubhshed in Ref. |27. 



The rest of the paper is organised as follows. In Section [TTl we present the details of our 
method, in Section IIIII our results and our conclusions are finally given in Section llVi 

II. METHOD 

In studies of static properties for diluted magnetic systems, a two step procedure is often 
used, combining density functional theory (DFT) and Monte Carlo (MC) techniques. We 
have here also used a two step procedure, but instead of a MC simulation as the second 
step, the time evolution of the spin dynamics has been explicitly simulated by treating the 
equation of motions for the atomic magnetic moments. We use atomistic spin dynamics 



(ASD) as derived in Refs. 



28 
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The inclusion of temperature in the simulations is treated 



with Langevin dynamics.— A thorough description of the details of our method is given in 
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Ref. 
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The central entity of the technique is the microscopic equations of motion for the 



atomic moments, mj, in an effective field, Bj, which is expressed as follows, 

= -7m, X [Bi + hi{t)] - 7— mi X (mi x [B^ + hi{t)]). (1) 

In this expression 7 is the electron gyromagnetic ratio and bj(t) is a stochastic magnetic 
field with a Gaussian distribution. The magnitude of that field is related to the damping 
parameter, a, in order for the system to eventually reach thermal equilibrium. We have 
used one heat bath and the time step for solving the differential equations was At = 0.01 
femtoseconds. 

The effective field, B,, on an atom z, is calculated from 

B. = (2) 

where only the part of the Hamiltonian, which represents interatomic exchange interac- 
tions, =^x, is considered in the present simulations. For this we use the classical Heisenberg 
Hamiltonian, 

where i and j are atomic indices and Jij is the strength of the exchange interaction, which 
is calculated from first principles theory. The exchange parameters Jij used h ere were 
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calculated by Kudrnovsky et al. from first principles theory as described in Ref. 
exchange interactions are strongly direction dependent and have been calculated up to and 
including the 39:th shell, which were all included in the present spin dynamics simulation 
(see Fig. [1]). Note in Fig. [1] that we have plotted the exchange parameters as a function 
of distance, not angle. As shown in Ref. 14| there is a strong angular anisotropy of these 
exchange interactions. The calculated Mn projected magnetic moment depends on the As 
antisite concentration,— and we have used the value of Ref. [ij, which is ~ /i = 4 /i^/atom. 



In principle the Mn projected moment is a little larger than 4, weakly increasing with As 
antisite concentration. Hence we consider As antisites to infiuence primarily the exchange 
interactions, not the size of the magnetic moments, a fact which is consistent with e.g. 
Refs. 



12 



and 



141 . Temperature effects of the size of the magnetic moment and strength of 



the exchange interactions, via changes of the electronic structure, have not been included in 
this work. The pair correlation function for equal times is defined as 
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Distance between atoms in units of lattice constants, r„/a 



FIG. 1: (Color online) Calculated exchange parameters for 5 % Mn doped in GaAs, over 39 atomic 
shells of the Ga sublattice, with y = 0.00 — 2.00 % the concentration of As antisites. The 39 shells 
have 29 unique distances, e.g. (0.5 0.5 2.0) and (0.0 1.5 1.5) both have distance ~ 2.12a. Data 



redrawn after Ref. 



14. 



Gij{t) = {uiiit) ■ m,(t)) 

= (m(ri,t) ■ m{r j,t)) 



(4) 



where (...) represents an average over all atoms in the simulation box. This quantity will be 
utilized to illustrate the dynamics as the system approaches thermal equilibrium. 



III. RESULTS 



A. Simulation of the Magnetisation 

To establish that a sufficiently large supercell was used, calculations were made for sys- 
tems with size L x L x L conventional fee unit cells, with L = 30, 40, 50. Periodic boundary 
conditions were imposed. With a Mn doping concentration of 5 %, this corresponds to 
5400, 12800 or 25000 magnetic atoms, respectively. The magnetization versus temperature 
curve, where the magnetization is the average of the magnetic moments in the supercell, is 
displayed in Fig. [2l It is noteworthy that convergence in the magnetization curve is achieved 
at L = 40 and for this reason all subsequent calculations were made with this or a higher 
value of L. We also note that Fig. [2] suggests an ordering temperature in the vicinity of 160 
K, which is close to the experimental value.- 

As the concentration of As antisites is increased, the interatomic exchange couplings 



TABLE I: Critical temperatures, Tc (K), for different As antisite concentration s, ?/ , given in per- 



cent. The Mn concentration is 5 %. MC stand for Monte Carlo data of Ref. 
current values from spin dynamic simulations. 



16 and ASD the 



As concentration 


MC 


ASD 


y=0.00% 


137 


160 


y=0.25% 




160 


y=0.75% 




135 


y=1.00% 


92 


115 


y=1.25% 


55 


75 


y=1.50% 


26 


38 



change from being purely ferromagnetic to having also an antiferromagnetic component. 
This has been shown to reduce the ordering temperature, both in experiment^, and in 
theory.— li^ii^ The temperature dependence of the magnetization from the ASD simulations 
is plotted in Fig. [31 for different values of the As antisite concentration y. We note that the 
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is the same. For a given set 



general behavior of the data in Fig. [3] and that of Refs. 
of pair exchange parameters ASD and MC simulations should indeed give the same result 
for the M(T) function and Tq. This follows as statistical mechanical considerations for 
thermal equilibrium require the energy of the magnetic degrees of freedom to be distributed 
according to a Boltzmann function^ in both methods. In Table. [T] our estimated values 
of the ordering temperature are seen to match fairly well with those calculated in Ref. [isl . 
The MC values for Tc were calculated using a fourth order cumulant method. Our ASD 
values for Tc were estimated from the points where the M{T) curves approaches zero. That 
circumstance might explain that our values are slightly higher than those in Ref. Q. 



B. Dynamical correlation 

The advantage of the spin dynamics approach over Monte Carlo simulations is that the 
true dynamics of the magnetism may be captured. As an example of this we show the 
time evolution of the macroscopic magnetic moment of our simulation box as a function of 
temperature, starting from a ferromagnetic (FM) or a completely random disordered local 
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FIG. 2: (Color online) Magnetization (normalized to Mq where Mq is the saturation magnetization 
at T = K), M{T), versus temperature, T, for Mn^jGai^^As with no As antisites, y = 0.00 %. 
L = 30, 40, 50 corresponds to different sizes of the supercell. 
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FIG. 3: (Color online) Magnetization M{T) (normalized to the saturation magnetization Mq at 
T = K) for Mn^^Gai-ajAs with As antisites. The antisite concentration ranges from y = 0.00 % 
toy = 2.00 %. 

spin configuration (rDLM). In the latter configuration, each spin is with uniform probability 
assigned a direction on the unit sphere, independently of the other spins. In Fig. |4]we show 
data for an As antisite concentration of y = 0.25 % at temperature T = 100 K. It should 
also be noted that the relaxation behavior is to some extent dependent on the damping 
parameter. In most of our simulations presented here we have used the value a = 0.10, but 
in Fig. Hlwe also show data for a = 0.03. Although in general smaller differences between 
results from the two damping parameters can be detected, we note that the main conclusions 
of our study are not influenced by this choice. A value of a = 0.10 is slightly higher than 
recent values for the damping in Mn^^Gai-^jAs, primarily obtained from experiments^^ but 
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also from theoretical estimates.— 

Figure m shows that the ferromagnetic starting configuration relaxes faster to the equilib- 
rium value, which is natural since the equilibrium configuration is closer to the ferromagnetic 
configuration than the rDLM configuration. The relaxation time is also different since the 
mechanisms to reach equilibrium are different. Starting from the ferromagnetic configuration 
the approach to equilibrium is primarily controlled by spin precession (magnon excitations), 
while starting from the rDLM configuration the approach to equilibrium involves both spin 
precession and spin flip processes. Since the relaxation in the two cases is different it is nat- 
ural that the relaxation rates are different. Figure H] also shows that the typical time-scale 
for the relaxation is ~ 5 — 10 ps for the ferromagnetic configuration and ~ 100 ps for the 
rDLM configuration at a T /Tc ratio 100/160, using a damping a = 0.1 and cell size L = 40. 
Graphs are also shown for a = 0.03 and it may be seen that in this case the time to reach 
equilibrium is approximately 250-300 ps. Hence when starting from a rDLM configuration 
the time to reach equilibrium seems to scale approximately linearly with the value of the 
damping parameter. 

A comparison can be made concerning the typical relaxation times for the system of 
interest here, Mn doped GaAs, and a typical ferromagnet, bcc Fe. When doing such a com- 
parison it is important to have similar values of the simulation cell, the damping parameter 
and the T/Tc, for the two systems, and we have indeed done this. The comparison was 
made at T/Tc= 0.62 (data for Fe not shown) and shows that the dynamics of Mn doped 
GaAs is slower with approximately 70 % than that of bcc Fe. To some degree this may be 
explained by the weaker local exchange field (B, in Eq. 1) of the Mn doped system, due to 
that it is a diluted system. 

The pair correlation function, Gij{t), depends not only on the distance between Mn 
atoms but also, as an immediate consequence of the directional dependence of the exchange 
interaction, on the direction in the lattice in which Gij{t) is calculated. In Figs. [SJlZlthe pair 
correlation function is plotted at different times as a function of distance between atoms. 
In these figures the deviation from a smooth decreasing behavior reflects indeed the angular 
anisotropy of the exchange interaction. Note that in these simulations we have used a 
temperature of 100 K and an As antisite concentration of 0.25 %, 0.75 % and 1.25 %. For 
the As 1.25 % antisite concentration the ordering temperature is below 100 K, but for the 
other As antisite concentrations the ordering temperature is above 100 K. The simulations 
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FIG. 4: (Color online) Time evolution of the average, normalized magnetization starting from 
ferromagnetic (blue) respective random (red) spin configurations for L = 40, with As Antisite 
concentration y = 0.25 % at temperature T = 100 K and with a damping parameter of 0.03. 
Similar simulations but with a damping parameter of 0.1 are shown in green and black for the 
random and ferromagnetic configurations, respectively. 
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FIG. 5: (Color online) Time evolution of the pair correlation function Gij{t) starting from ferro- 
magnetic (blue) respective random (red) configurations and with L = 40 in the simulation. Values 
are obtained for y = 0.25 % and T = 100 K. The equilibrium pair correlation is shown in black. 

were started both from a rDLM and a FM configuration. A few general remarks can be 
made concerning the data in Figs. 013 For instance, nearest neighbor (NN) spins reach their 
equilibrium values faster (typically after 1 ps), whereas with increasing distance the dynamics 
of the correlation slows down (5-10 ps). It is also clear that for Figs.[5]and [6]the FM starting 
configuration seems to reach equilibrium faster than the rDLM configuration, whereas for 
Fig. [7] the rDLM configuration reaches equilibrium faster than the FM configuration. It 
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Distance between atoms in units of iattice constants, r./a 



FIG. 6: Same as Fig. [5] but for y = 0.75 




Distance between atoms in units of lattice constants, r,./a 



FIG. 7: Same as Fig. [5] but for y = 1.25 %. 



is natural that the systems shown in the former two figures reach equihbrium from the 
ferromagnetic configuration faster, since these systems are below the ordering temperature, 
where a ferromagnetic starting point should be closer to the equilibrium configuration. We 
also note that for all systems there is a very strong short ranged order, involving several 
neighboring atomic shells, which is present at all temperatures, even up to and above the 
ordering temperature (cf. in Fig. [7]). 

Turning to the details of the curves we note that e.g. in Figs. [5] and [6] the NN relaxation 
time is close to 1 ps, both for the rDLM and FM starting configurations. The long ranged 
(39 shell distance) relaxation is ~5 ps for the FM starting configuration and > 10 ps for 
the rDLM configuration. For the system which is above the ordering temperature a slightly 
different behavior is observed. In Fig. [7] the rDLM configuration reach equilibrium after 1 
ps, for the NN and longer ranged distances. The FM configuration is here relaxing slower. 
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requiring 5 ps. To some degree, the dynamics above the ordering temperature is similar to 
that below the ordering temperature. 

The data in Figs. [MZl show that as the As antisite concentration increases, antiferromag- 
netic, superexchange interaction between the Mn- atoms becomes more dominant. This is 
especially obvious for the NN interaction. This result is consistent with the analysis pre- 
sented in Ref. |3J where a non-collinear coupling between NN Mn atoms was found when 



competing ferromagnetic and antiferromagnetic interactions play role. We have plotted the 
distribution P{6) of angles between the local moments and the average magnetization in the 
upper panel of Fig. Mioi y = 1.25% and y = 1.50% at T = K. The distribution of angles 
sums up to an average magnetization of M/Mq = 0.87 for y = 1.25% and M/Mq = 0.61 
for y = 1.50% (cf. Fig. [3]). To analyze P{0) the distribution of angles for the disjoint sets 
of moments with no NN:s respective at least one NN are plotted. In the middle panel of 
Fig. [8] only the moments with at least one NN have been considered. The antiferromagnetic 
coupling between NN:s makes sure that a non-negliable fraction of these moments are close 
to being antiparallel with the direction of the average magnetization. In the lower panel are 
plotted P{6) for local moments that do not have any NN:s. With no strong antiferromag- 
netic coupling acting on these moments their angles towards the average magnetization is 
effectively confined to the intervals - 25° {y = 1.25%) and - 45° {y = 1.50%). 

Our distributions P{0) for moments with at least one NN can be compared with Fig. 5 



in Ref. 



34j . For antisite concentration y = 1.25% the ASD simulations have a peak in the 



intervall of 4 — 6° where the canted spin model (GSM) instead has a delta peak for 6 = 0°. 
The mean angle of the respective distribution differ substantially {Oasd = 41°, 9csai ~ 30° 
). For antisite concentration y = 1.50% the difference in mean angles is small {Oasd ~ 82°, 
0c SM ~ 83° ) although the shape of the distributions are different. The very same exchange 



34. We 



couplingsi^ were used for the present ASD simulations and the CSM study in Ref. 
conclude that the various approximations in the two schemes account for the differences in 
the histograms of the angles. 

In addition, the possibility of non-collinear magnetism in zinc-blende MnAs was discussed 
in Ref. 



35l . which may be important for understanding samples with an inhomogeneous 
distribution of Mn atoms where Mn rich regions are created. Inside such a region, where 
the composition locally is close to MnAs, non-collinear magnetic couplings between the Mn 
atoms are likely to occur. 
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Canting angle (degrees) 

FIG. 8: (Color online) Canting for y = 1.25% (green) respective y = 1.50% (blue) at K. The 
upper panel shows the distribution of angles for all moments in the simulation cell. The middle 
panel shows the distribution for moments with at least one NN. Data for the disjoint set of moments 
with no NN:s are shown in the lower panel. 

The distribution of the individual scalar products mj ■ mj (which we will refer to as 
the spin-product) for a simulation at 100 K and with a defect concentration of 0.25 % is 
shown as histograms in Figs. [9|T0] and fTTl for times t = 1,5 and 10 ps, respectively, from a 
simulation which starts from a completely random configuration of spins. The histograms 
were calculated as follows. At a given time, mj ■ nij was calculated between a given selected 
Mn atom and all the Mn atoms in the nearest neighboring shell (shell 1), and between the 
chosen Mn and all Mn atoms in the next nearest neighboring shell (shell 2), etc. This was 
then repeated over all Mn atoms in the simulation box. In this procedure all spin-products 
are summed twice and we have corrected for this fact. The number of spin-products with a 
given value are then shown in histograms Figs. [QlfTTl for various times in the simulation. 

Figures [QlfTTl show that after 1 ps the NN spin-product is equilibrated to its saturation, 
whereas spin-products between Mn atoms further apart (shell 2-10) have not quite reached 
equilibrium, but is nevertheless deviating strongly from completely disordered spins (which 
would show a uniform distribution). An almost uniform distribution is the general behavior 
for shell 39. After 5 ps into the simulation (Fig. [TOl) the spin product between almost all 
shells have reached their equilibrium distribution and deviate very little from the data shown 
for 10 ps (Fig. [TTj) . The general behavior of Figs. [QlfTTl is rather similar to a droplet model 
of magnetic correlations, e.g. as used to analyze spin-glass materials, in that a correlation 
around a given spin grows outwards in time after the simulation has starteGp^*^"^. 
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FIG. 9: (Color online) Histogram of the spin product function lUi ■ nij starting from random 
configurations after 1 ps. The simulation box had L = 40. The temperature is 100 K and the 
As antisite concentration is 0.25 %. The value of the spin product is given on the x-axis and the 
number of atoms with a given spin product is given on the y-axis. The pair correlation between 
different atomic shells is marked in each panel with the shell number in question (e.g. shell 1, shell 
2 etc.). The grading of the y-axis varies with the number of atoms in the respective shell. 
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FIG. 10: (Color online) Same as Fig. [9] but after 5 ps. 

The data in Figs. ISlfTTl suggest that local relaxations cease after some 10-20 ps for a rDLM 
configuration. These relaxation times are some ten times faster than what Fig. H] suggests for 
the same configuration. The reason for this difference is that the data in Figs. [5lfTT] measure 
time scales of local correlations, whereas the data in Fig. H] reflects the whole simulation cell. 
Taken together, the data in Figs. l4lfTT] suggest that local spin-droplets develop after some 
10-20 ps, and that the correlation between such droplets takes 5-10 times longer to develop. 
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FIG. 11: (Color online) Same as Fig. [9] but after 10 ps. 



C. spin-glass analysis, aging and autocorrelation 



The possibility of spin-glass behavior among DMS materials, as discussed in the Intro- 
duction, naturally calls for a theoretical analysis, and the spin-dynamics simulations used in 
this work is a good method for undertaking such an analysis. Hence, from our simulations 
we have studied aging phenomena, which are typical for spin-glasses, on 5 % Mn doped 
GaAs with various concentrations of As antisites. The reason for using the As antisite con- 
centration as a parameter is that one gradually increases the amount of anti-ferromagnetic 
interactions in the lattice, which starts from purely ferromagnetic interactions, when the As 
antisite concentrations is zero (see Fig. [1]). 

A slowly relaxing system has physical observables that break time translation invariance 
before it eventually equilibrates. Two-time observables can sample the characteristics of the 
approach to equilibrium. The spin autocorrelation function 

Co{t^ + t, t^) = (mi(t^) • nii{t^ + t)) (5) 

is the simplest two-time quantity that has been used to study the aging of spin-glasses^^. 
The autocorrelation is not directly measurable in experiments, where instead the closely 
related zero (in practice low) field dynamic susceptibility Wli{tyj,t)/h has been used^. The 
calculation of Co{tw + t^tw)-, was done by starting a simulation from a completely random 
configuration of Mn moments. A waiting time, t^;, was selected and the scalar product 
between a given moment mj(t^) at this time and the same moment a time t later (i.e. 
^i{^w + ^) ) was evaluated. Note that Eq. [5] also contains an average over all magnetic 
atoms of the simulation box. 
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FIG. 12: (Color online) Autocorrelation Co{t^ + t,tyj) for 5 % Mn doped GaAs, with 0.25 % As 
antisite concentration. Simulations were made at T = 10 K, and a = 0.03 (top panel) and a = 0.10 
(bottom panel). 

Before entering the details and the analysis of the autocorrelation functions of Mn doped 
GaAs, shown in Fig. [12] (for 0.25 % As antisites) , Fig. [13] (for 1.75 % As antisites) and 
Fig. [14] (for 2.00 % As antisites), we note that data shown for Mn doped GaAs are obtained 
from averages over up to 20 simulations. In making these averages we used 4 different 
configurations of Mn atoms, and we used 5 different seeds for generating the heat baths used 
in the simulations. It was found that making an average over different configurations was 
more important for obtaining reliable curves to analyse. For comparison, the autocorrelation 
of a ferromagnetic system (bcc Fe) below and above the ordering temperature is discussed 
briefly (Fig. [T5|) . 

The choice of concentrations of As antisites was made to investigate the dynamics as 
a function of increasing antiferromagnetic interactions, in samples which have a net ferro- 
magnetic moment (Figs. [12] and [T3|) and to compare this behavior with a system for which 
there is no long ranged order (but with finite atomic Mn moments. Fig. [T4l) . The data in 
Figs. [T2][T3] and [H] show that the main characteristics of the autocorrelation functions stud- 
ied here are not influenced by the choice of damping parameter. The figures also show that 
for the systems with finite magnetization at T = 10 K (Figs. [T2] and [T3|) the autocorrelation 
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FIG. 13: (Color online) Same as in Fig. [12] but for 1.75% As antisites. 

function approaches a finite value witli increasing time, t, for long enough waiting times, tw 
The data in Figs. [12] and [13] should be compared to a normal ferromagnetic material, like bcc 
Fe, and for this reason we show autocorrelation functions for this well-known ferromagnet at 
T=100 K in Fig. [15] (upper panel). Note that for bcc Fe all autocorrelation functions reach 
a saturation value faster than it does for Mn doped GaAs, but that otherwise the shape of 
the autocorrelation functions are not too different. 

In section IIII Bl we discussed how the system relaxes from a rDLM or FM configuration 
to reach a thermally equilibrated configuration. This discussion was based on the values 
of the spatial correlation function at different times and will here be extended with the 
features that can be seen in the autocorrelation curves. For long enough waiting times 
(curve marked equil in the lower panel of Fig. [12] and tw > 10 ps in Fig. [T5|) . the system has 
already reached thermal equilibrium. Here the autocorrelation, starting from unity at an 
observation time t = 0, quickly settles to a value that equals the value of the magnetization 
the system posseses at that temperature. The individual spins fluctuate and precess around 
the direction of the effective magnetic field and their individual autocorrelation functions 
fluctuate correspondingly. Averaged over the ensemble these fluctuations vanishes and the 
autocorrelation function, after the initial decay from unity, is constant. For short waiting 
times {0 < tyj < 2.5 ps in the lower panel of Fig. [12] and < t^o < 0.625 ps in Fig. [T5l) . the 
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autocorrelation quickly drops from unity to essentially zero. 

The more interesting features are to be found for intermediate waiting times (5 < < 
200 ps in the lower panel of Fig. [Inland 1.25 < t«, < 5 ps in Fig. [T^ . The autocorrelation 
initially decays to a value equal or slightly smaller than the value for the longer waiting 
times. With increasing observation time the autocorrelation eventually drops from this level 
and starts a damped oscillation. The value around which the autocorrelation oscillates cor- 
responds to the value the magnetization has reached for a specific waiting time. When 
comparing the frequencies of the oscillations in Figs. [12] and [15] (which represent an average 
over all spins in the simulation cell) with the typical frequencies of individual spin trajec- 
tories, we observe that they are of comparable magnitude. The observed oscillations can 
therefore be explained from the fact that each atomic spin precesses in its local exchange field 
and that this spin precession is faster than the rate at which the local exchange field reaches 
its equilibrium orientation. Summing over all spins in the simulation box then gives rise to 
the oscillatory behaviour in Figs. [12] and [15] The damping term in Eq. 1 ensures that the 
oscillations are damped. If one would make an average over different configurations of the 
Mn atoms and different heatbaths, the oscillation diminish in magnitude with the number 
of simulations. This is illustrated in the middle panel of Fig. [15], where the oscillations are 
seen to have been strongly suppresed. 

For Mn doped GaAs with 2.00 % As antisites (Fig. UM the situation is very different 
compared to the simulations with lower As antisite concentration, since the autocorrelation 
function approaches zero, for all waiting times. The reason is that the sample for T = 10 K 
is paramagnetic. We note however that there are signs of slow dynamics in the simulations, 
since the different autocorrelation functions do not lie on top of one another, but are spread 
out over a time interval which covers approximately one decade in time. This is distinctly 
different from bcc Fe in the paramagnetic phase, at temperatures well above the ordering 
temperature (see Fig. [15]- lower panel), where all autocorrelation functions lie on top of one 
another, independent on waiting time. 

For very low temperatures, one can expect a phase transition also for Mn doped GaAs 
with 2.00 % As antisites. For a pure spin-glass system, the spin dynamics will exhibit critical 
slowing down on approaching the transition temperature from above. Below the spin-glass 
transition temperature, the different autocorrelation functions should not approach zero for 
any time t, and curves corresponding to different waiting times should never lie on top of 
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one-another for any extended interval of time. Our simulations for T = 2 K (data not 
shown) do indeed show a very slow dynamics. For instance, starting from a ferromagnetic 
configuration or from a rDLM configuration we do not reach equilibrium within the time of 
the simulation. However, for computer simulations made with limited cell size, the system 
will eventually exhaust the length scales allowed by the periodic boundary conditions and 
the autocorrelation function will eventually go to zero. In Ref. |39| Berthier et al. treats the 
mechanism of finite size effects in Monte Carlo simulations on a archetypical Heisenberg spin- 
glass. Problems connected to finite size are also present in our spin dynamics simulations, 
and we are left to conclude that our data on 2.00 % at T = 2 K are not sufficient to allow 
for any definitive statements on the nature of the phase transition. 

It is meaningful to compare the results of Figs. [T^fT^ and [HI with simulations of typical, 
well known spin-glass systems, like Mn doped Cu (fee). We have undertaken simulations 
on Mn doped Cu (data not shown) and we find that these simulations do indeed show 
typical spin-glass characteristics, but details of this simulation are outside the scope of the 
present analysis and will not be discussed further here. Nevertheless, we conclude that our 
simulation package results in the expected behavior for well established spin-glass systems. 
The fact that Mn doped GaAs, for several concentrations of As antisites, with varying degree 
of competing ferromagnetic and antiferromagnetic exchange interactions, does not display 
spin-glass characteristics suggest that it is incorrect to label this system as a spin-glass. 



D. Visualization 



Scientific visualization techniques enables further insight into spin dynamic simulations. 
We have produced a set of movies that are accessible on the webpage of Ref. These 
simulations were done without As antisites at a temperature of 100 K (200 K) which is 
below (above) the ordering temperature. The movies reveal that also above the ordering 
temperature the spins are correlated on short distances. The difference occurs for longer 
distances, where for 100 K, the spins are still correlated, but for 200 K they are not. 
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FIG. 14: (Color online) Same as in Fig. [I2] but for 2.00% As antisites. 
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FIG. 15: (Color online) Autocorrelation Co{tw + t,tw) for bcc Fe at T = 100 K for a single run 
(upper panel) and average over 20 runs (middle panel). For comparison also the autocorrelation 
at r = 1500 K is shown (lower panel). In the simulations we had a = 0.10. 
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IV. CONCLUSIONS 



In conclusion we have applied atomistic spin dynamics simulations to 5 % Mn doped in 
GaAs, with and without As antisite defects. It should be noted that also other types of 
defects have been found experimentally for this system, e.g. Mn interstitials. Both types 
of defects introduce antiferromagnetic interactions in the lattice and for simplicity we have 
here only considered the As antisites. Our results show ordering temperatures and a general 
behavior of the magnetization curve which are in agreement with Monte Carlo simulations 
and experiment. 

The dynamical response has also been investigated and we find that the dynamics between 
NN Mn atoms is, as revealed by the pair correlation function, considerably faster than 
the more long ranged interaction. In addition we find that short ranged order exists up 
to and even above the ordering temperature. We also find that starting from random 
configurations, the pair correlation function, for distances up to 4a (a = lattice parameter), 
reaches equilibrium on a time scale of 10-20 ps. The dynamical response, as revealed by the 
autocorrelation function, shows that Mn doped GaAs, does not display spin-glass behavior, 
for moderate {y = 0% — 1.75%) concentrations of As antisites. For the highest defect 
concentration, y = 2.00%, the nature of the phase transition (ferromagnetic or spin-glass) 
of the simulation cell cannot be accurately determined from our data. 

The present article focuses on one III-V semiconducting host, GaAs, and one magnetic 
element, Mn, in the low concentration regime. However, the general behavior of the mag- 
netism of most doped semiconductors is basically similar, with an exchange interaction 
which decays exponentially with distance due to the presence of a band gap in the system, 
overlapped with an angular anisotropy of the interatomic exchange, and the possibility to 
introduce competing ferromagnetic and antiferromagnetic interactions via the defect con- 
centration. For this reason we believe that the main results concerning the dynamics of 
the presently studied system should be representative for several DMS materials, at least 
in the diluted limit of magnetic dopants, and that they should be characterized as slow 
magnets. A notable exception from this conclusion may be the II- VI semiconductors, where 
the solubility of magnetic atoms can be considerably higher (up to ~ 30 %) than what 
we investigate here. For some of these materials2ii22i2^ experimental data do indeed show 
spin-glass behavior. Investigations of such systems are underway. 
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